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The electromagnetic charge operator in a two-nucleon system is derived in chiral effective field 
theory (^EFT) up to order eQ (or N4LO), where Q denotes the low-momentum scale and e is 
the electric charge. The specific form of the N3LO and N4LO corrections from, respectively, one- 
pion-exchange and two-pion-exchange depends on the off-the-energy-shell prescriptions adopted for 
the non-static terms in the corresponding potentials. We show that different prescriptions lead to 
unitarily equivalent potentials and accompanying charge operators. Thus, provided a consistent 
set is adopted, predictions for physical observables will remain unaffected by the non-uniqueness 
associated with these off-the-energy-shell effects. 



PACS numbers: 12.39.Fe, 13.40.-f 
I. INTRODUCTION 

The chiral symmetry exhibited by quantum chromody- 
namics (QCD) dictates that the pion couples to baryons, 
such as nucleons and A-isobars, by powers of its momen- 
tum Q. As a consequence, the Lagrangian describing 
these interactions can be expanded in powers of Q/A x , 
where A x ~ 1 GeV specifies the chiral-symmetry break- 
ing scale. Thus, classes of Lagrangians emerge, each char- 
acterized by a given power of Q/A x and each involving 
a certain number of unknown coefficients, so called low- 
energy constants, which are then determined by fits to 
experimental data (see, for example, the review papers [l[ 
and and references therein). This approach, known 
as chiral effective field theory (%EFT) , can be justifiably 
argued to have put low-energy nuclear physics on a more 
fundamental basis by providing, on the one hand, a direct 
connection between the symmetries of QCD — in partic- 
ular, chiral symmetry — and the strong and electroweak 
interactions in nuclei, and, on the other hand, a practical 
calculational scheme susceptible, in principle, of system- 
atic improvement. 

The model for the nuclear electromagnetic current — 
the space-like part of the four current — in xEFT up to 
one loop was derived originally by Park et al. Q, us- 
ing covariant perturbation theory. In the last couple 
of years, two independent derivations, based on time- 
ordered perturbation theory (TOPT), have appeared in 
the literature, one by the present authors [J] and the 
other by Rolling et al. [!|. There are technical differ- 
ences in the implementation of TOPT, which relate to 
the treatment of reducible diagrams and are documented 
in considerable detail in the above papers. However, the 
resulting expressions in Rcfs. Q and [f| for the two-pion- 
exchange currents (the only ones considered by the au- 



thors of Ref. are in agreement with each other, but 
differ from those of Ref. [3[ , in particular in the isospin 
structure of the Ml operator associated with the one- 
loop corrections — see Pastore et al. (2009) Q for a com- 
parison and analysis of these differences. 

Rolling et al. also provided the first treatment of loop 
corrections to the nuclear charge operator — the time-like 
part of the four current — associated with two-pion ex- 
change mechanisms. Of course, there had been earlier 
studies of the two-nucleon charge operator, notably those 
of Refs. [gUl], but they had been limited to its isoscalar 
component, and therefore had only retained tree-level 
corrections (the two-pion-exchange loop contributions 
are isovector) . These earlier studies also included predic- 
tions for the charge and quadrupole form factors of the 
deutcron, which were in reasonable agreement with data 
obtained from measurements of elastic electron-deuteron 
scattering cross sections at low-momentum transfer. 

The primary objective of the present work is to extend 
the formalism developed in Rcfs. [H, to derive system- 
atically the two-nucleon electromagnetic charge opera- 
tor in vEFT, including up to one loop corrections. As 
we shall see below, this is not a straight-forward task, 
since the derivation of such an operator necessarily entails 
the study of non-static contributions to the one-pion- 
exchange (OPE) and two-pion-exchange (TPE) poten- 
tials. In the OPE sector, this inter-connection between 
non-static contributions and the charge operator was in- 
vestigated long ago by Friar [l(| in the context of a Foldy- 
Wouthuysen reduction procedure and a time-dependent 
perturbation theory, which consistently retained correc- 
tions up to order (v/c) 2 . In particular, Friar showed that 
i) the charge operators so derived depend on the specific, 
but arbitrary, off-the-energy extension — i.e., on the cor- 
rections beyond the static limit, such as those induced by 
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retardation effects adopted for the OPE potential, and 
that ii) these different operators (and corresponding OPE 
potentials) are related to each other by a unitary trans- 
formation and, therefore, their intrinsic lack of unique- 
ness has no consequence on the predictions for physical 
observables. 

In this paper, we examine these issues from a %EFT 
perspective. We are interested in constructing the charge 
operators up to one loop and, hence, need to include 
non-static corrections not only in the OPE, but also in 
the TPE potential. We show that the resulting oper- 
ators, while not unique because of the off-thc-encrgy- 
shell ambiguity referred to above, are nevertheless re- 
lated to each other by a unitary transformation. Thus 
the present study puts Friar's original considerations in 
the modern framework of %EFT and extends them to 
the TPE sector. In Sec. HD and App. H we list those 
terms in the chiral Lagrangians (and corresponding in- 
teraction Hamiltonians) that are relevant to our purpose 
here. In Scc. lIIII wc provide an overview of the derivation 
of the OPE and TPE potentials and charge operators. 
Specifically, we explore the connection between the am- 
plitude calculated in xEFT and the strong and electro- 
magnetic potentials, which are derived from it and are 
used in quantum-mechanical formulations, based on the 
Lippmann-Schwinger or Schrodingcr equations. Power 
counting allows us to establish a criterion to make this 
connection precise. Elsewhere [1, @, we have referred 
to the latter as "accounting for recoil-corrected reducible 
contributions." The present formulation is especially apt 
to shed light on the inter-dependence between charge op- 
erators and potentials (at the OPE and TPE level), and 
the associated ambiguity arising from off-thc-energy-shcll 
extrapolations prescribed for the latter. 

In Sec. IIIII we also provide an explicit expression for 
the unitary transformation, and show, in particular, that 
different (non-static) versions of the TPE potential are 
unitarily equivalent. As mentioned above, in the OPE 
sector this result has been known for a while [HI E| ■ 

Section IIVI contains a summary of the derivation of 
the two-nucleon charge operators up to order e Q in- 
cluded, or next-to-next-to-next-to-leading order (N3LO 
in short), where e is the proton electric charge and Q 
denotes generically the low-momentum scale. No loops 
enter at this order. The different forms of the OPE charge 
operator at N3LO exhibit the same unitary equivalence 
as the non-static corrections to the OPE potential at or- 
der Q 2 (N2LO)— this too is well known [Io|, EH, albcit in 
a different context. 

In Sec. [V] we discuss the static one-loop corrections to 
the charge operator at N4LO (eQ). In particular, we 
list those corresponding to two different off-the-energy- 
shell prescriptions for the OPE and TPE potentials, and 
show that they are unitarily equivalent. A fairly detailed 
account of their derivation is provided in App. [Bj The 
loop integrals entering the individual terms at N4LO are 
ultraviolet divergent; however, their sum is finite, in par- 
ticular it vanishes in the limit in which the momentum q 



carried by the electromagnetic field is zero. This was 
to be expected, since i) symmetry arguments prevent 
the presence of counterterms at this order (and there- 
fore the possibility of re-absorbing ensuing divergencies 
into them), and ii) charge conservation demands that at 
q = the charge operator merely counts the number of 
charged particles (i.e., protons) in the system — a require- 
ment already fulfilled at LO. Summary and conclusions 
are presented in Sec. lVIl while, for future convenience, we 
give the configuration space representation of the N4LO 
charge operators in App. [Cl 



II. RELEVANT INTERACTION 
HAMILTONIANS 

Here we only list the interaction Hamiltonians relevant 
for the derivation of the nuclear electromagnetic charge 
operator up to order e Q (N4LO), see App. [Al for notation 
and a summary of the corresponding Lagrangians (l2l . 

H-kn = /dxA+ [2y±<r ■ V^a + ^ • (tt x 9°7r)]7V,(l) 
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The resulting vertices behave, relative to the low- 
momentum scale Q, in the following way: H^n ~ Q\ first 
term in H 7 n ~ e Q °, remaining ones ^ eQ 2 ; H 7n ^ eQ; 
first term in H^^n ~ eQ, second one ~ eQ 2 . There is 
also a contact interaction, 

H C t = - I 



C s (N^N) (N^N) 



+ C T (N^ctN) ■ (N^erN) 



(5) 



which enters the derivation of the N4LO charge operator. 
The accompanying vertex scales as Q . 



III. FROM AMPLITUDES TO POTENTIALS 

We begin by considering the conventional perturbative 
expansion for the two-nucleon (NN) scattering ampli- 
tude 



(f\T\i) = (f\H 1 
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where | i) and | /) represent the initial and final two- 
nucleon states of energy Ei — Ef, Ho is the Hamilto- 
nian describing free pions and nucleons, and Hi is the 
Hamiltonian describing interactions among these parti- 
cles (Sec.|TT|. The evaluation of this amplitude is carried 
out in practice by inserting complete sets of Hq eigen- 
states between successive terms of H±. Power count- 
ing is then used to organize the expansion in powers of 
Q/A x <C 1, where A x ~ 1 GcV is the typical hadronic 
mass scale. 

In the perturbative series, Eq. a generic (reducible 
or irreducible) contribution is characterized by a certain 
number, say N, of vertices, each scaling as Q ai x Q~^^ 2 
(i=l, . . . , N), where oti is the power counting implied by 
the relevant interaction Hamiltonian and /?,; is the num- 
ber of pions in and/or out of the vertex, a correspond- 
ing N—l number of energy denominators, and possibly 
L loops Jj|. Out of these N-l energy denominators, 
Nk of them will involve only nucleon kinetic energies, 
which scale as Q 2 , and the remaining N — Nk — 1 will 
involve, in addition, pion energies, which are of order Q. 
Loops, on the other hand, contribute a factor Q 3 each, 
since they imply integrations over intermediate three mo- 
menta. Hence the power counting associated with such 
a contribution is 



equation, 



f N 
\i=l 



-(N~N K -1) q-2N k 
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Clearly, each of the N — Nk — 1 energy denominators can 
be further expanded as 



Ei — Ej — oj^ 
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1 



Ei — Ej 



(Ei — Ei ) 
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where Ej denotes the kinetic energy of the intermediate 
two-nucleon state, and the pion energy (or energies, 
as the case may be) — the ratio (Ei — Ei)/u> n is of order 
Q. 

The Q-scaling of the interaction vertices and the con- 
siderations above show that T admits the following ex- 
pansion: 



rp _ Jl(0) , Jl(l) _|_ Jl(2) , 



(9) 



where ~ Q n . For example, the time-ordered di- 

agrams contributing to and are illustrated in 
Fig. [T] where the pion line (pion line with a crossed 
circle) indicates that only the leading — (next-to- 
leading — (Ei — Ej)/uj 2 ) term is retained in the expansion 
of the associated energy denominator, Eq. ©. Except for 
App.[Bl this notation will not be used any further below, 
but it is understood that energy denominators involving 
pions are expanded as in Eq. ([8]). 

Our objective is to derive a two-nucleon potential v 
which, when iterated in the Lippmann-Schwinger (LS) 



v + v Gq v + v Gq v Gq v + 



(10) 



leads to the T-matrix in Eq. ([5]), order by order in the 
power counting. In practice, this requirement can only be 
satisfied up to a given order n* , and the resulting poten- 
tial, when inserted into the LS (or Schrddinger) equation, 
will generate contributions of order n > n* , which do not 
match TW. In Eq. (|10p . Go denotes the free two-nucleon 
propagator, Go = — Ej + irf), and we assume that 



U = W (0) +U (D +U (2) 
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where the yet to be determined »W is of order Q n . We 
also note that, generally, a term like t/ m ) Go v 1 ^ 1 ' is of 
order Q m+n+1 , since Go is of order Q~ 2 and the implicit 
loop integration brings in a factor Q 3 . Having established 
the above power counting, we obtain 
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(15) 



where the "recoil-corrected" two-nucleon potential 

explicitly constructed in Refs. [3, [9| up to order n = 2, 
or N2LO. The LO term, u'"', consists of (static) one- 
pion-exchange (OPE) and contact interactions, while the 
NLO term, vanishes, since the contributions of di- 
agrams (d) and (e) in illustrated in Fig. [TJ add up 
to zero, while the remaining diagrams represent itera- 
tions of v(°\ whose contributions are exactly canceled by 
[v^ Gqv^] — complete or partial cancellations of this 
type persist at higher (n > 2) orders. The N2LO term, 
which follows from Eq. (|14p . contains two-pion-exchangc 
(TPE) and contact (involving two gradients of the nu- 
cleon fields) interactions. It is derived in Ref. How- 
ever, there is a recoil correction of order n = 2 to the 
OPE potential, which was ignored in that paper. In mo- 
mentum space, it is given by 



i4 2 >(i/ = 0) = «<°>(k) 



(E{-E 1 ) 2 + (EJ l -E 2 ) 2 
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where is the LO OPE potential, 

^ w 2iJlfiJE, 
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k = pi p'j = p' 2 — P2 is the momentum transfer, and 
Pj and Ej (p^- and Ej) are the initial (final) momentum 
and energy of nuclcon j. Obviously, on the energy shell 
Ei = Ef implicit in Eq. (fT4")h the above expression is 
equivalent to one in which, for example, 



l) = -«( )(k) 



{E{ - Ei) (E< - E 2 ) 



(18) 



(2) 

In fact, there is an infinite class of v\ jv) corrections — 
labeled by the parameter v [0, [Hj] — which, while 
equivalent on the energy-shell, arc different off the 
energy-shell, and therefore lead to different potentials 
v^(y) in Eq. (fT3|) . Indeed, for the choices of v^(v) 
in Eqs. (|T6|) and (fT8|) the corresponding corrections to 
the TPE term (from direct and crossed box diagrams), 

M> read 



<4tV = o) 



(3 + 2ti-t 2 )/ (<ti • q 2 )(eri • qi) 

J qi 



2F* 

x (er 2 ■ qi)(cr 2 • q 2 ) 



qi 

Ei — Ei 



E 2 - E' 2 
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and 

4°V = i)=4 3 > = o) + 

x (E 1 + E 2 



(19) 



1 / 4°)(q 2 )^(qi) 

Z Jqi 



2 

Ei - E 2 



,(20) 



where, as indicated in panel (a) of Fig. |4] (App. |Bj, qi 
and q 2 (u>i and w 2 ) are the momenta (energies) of the 
two exchanged pions (with q 2 = k — qi), Ej and Ej are 
the intermediate nucleon energies, and 



ds 



(21) 



Friar [H HI, and later Adam et al. [TO], have ar- 
gued that the different off-the-energy-shell extrapolations 
(v) are unitarily equivalent. (See also Ref. []1| for the 
implication of the unitary equivalence on the TPE three- 
nucleon potential.) We show below that, as a matter of 
fact, this unitary equivalence remains valid for v^J (v) 
(as well as for the electromagnetic charge operators at 
one loop), thus extending the results of the authors of 
Refs. [Ill El to the TPE sector. Up to order n = 3 (i.e., 
Q 3 ), the two- nucleon Hamiltonian in the center-of-mass 
(CM) frame can be written in momentum space as 

H{v; p',p) = K<--Q(p', p) + 4°)(p'- p) + w&V - p) 

+ ^ 2 W,p)+« 2 'W,p) , (22) 

limiting our considerations to OPE and (box) TPE po- 
tentials only. There are, of course, v ^' terms originating 
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FIG. 1: Time ordered diagrams illustrating the contributions 
to the T (0) , panels (a)-(c), and T (1) , panels (d)-(n), AW 
scattering amplitudes. Nucleons and pions are denoted by 
solid and dashed lines, respectively. Pion lines (pion lines 
with crossed circle) indicate that only the leading Q ~ (next- 
to-leading Q °) term in the expansion of energy denominators, 
Eq. |[8]), are retained in the corresponding amplitudes. See 
text for explanation. 



from higher order chiral Lagrangians fl2l fl3l | , but these 
have no relevance for the discussion to follow, and are 
therefore ignored below. In Eq. (|2"2"|) . K^ 1 ^ denotes the 
kinetic energy term of order n = — 1, 



K^(p',p) = (2- ! r) 3 6(p'-p)p 2 / 



TO AT 



(23) 



(2) 

mjv is the nucleon mass, and v\ (y) has been given in 
Refs. [HQ as 



v™{v;p',p) = {l-2v) 



^°V-p) [ P i2 -p 2 f 



(p' - p) 



Am 2 N 



(24) 

which the v = 0, 1 expressions listed above reduce to (in 
the CM frame). These Hamiltonians are related to each 
other via 



H(y) = e- iU{v) H(v = 0) c +lC/(l/) , 
where up to NLO the operator i U(v) is 

* U{y; p', p) ~ i (1/; p', p) + i U^{y; p', p) 
with 



i^°)( V ;p',p) = -i/ 7 ^ 



^ 0) (P'-P) P' 2 ~P 2 



(p' - p) 2 + ml 2 Wiv 



:uM(y; p',p) 



" /" «l 0) (p' - s)"x u; (s - p) 



(0), 



(25) 



(26) 



(27) 



(28) 



2 7 S (p' - s) 2 + m 2 
The unitary equivalence up to order n — 3 implies 

= H (v = 0) + [k^+v^, iU {0 \u) 

+ [ir^ 1 ), ii/w , (29) 
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since each commutator brings in an additional factor Q 3 
due to the implicit momentum integrations. A direct 
evaluation with v = 1 shows that H(v = 1) ensues, in- 

eluding i4 {y = 1) and \ v = -0 as gi vcn m Eq s - CHI) 
and (f^Oj) — note that in the CM frame qi = p s, 

q2 = s — p', E1+E2 = s 2 /m;y, and s is the loop momen- 
tum. Both these ^-dependent corrections are relevant 
for the derivation of the nuclear charge operator up to 
N4LO, to which we now turn our attention. 

The electromagnetic interactions are treated in first 
order in the perturbative expansion of Eq. ([6]), and the 
transition operator can be expanded as 



(30) 



where T 7 is of order e Q n (e is the electric charge) . The 
nuclear charge, p, and current, j, operators follow from 
v 7 = A p — A ■ j, where = (A° , A) is the electromag- 
netic vector field, and it is assumed that i> 7 has a similar 
expansion as T 7 . The requirement that, in the context 
of the LS equation, t> 7 matches T 7 order by order in the 
power counting implies the following relations: 



,(-1) 



T (-3) 
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T (-2) 
7 
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v ( ~ 2) G w (0) G w (0) + permutations 



vi~V Gov^ + vW Go^- 1 ) 



Gov^+v^ Go „(-3) 
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v\~V G vM G v^ GoU (o) 
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(34) 



+ permutat 



10ns 



vi~ 3) G v {2) G w (0) + permutat 



10ns 



(35) 



where v^ l) = A p^ - A-jW, arc the AA potentials 
constructed in Eqs. ((TS]) - (fT5]) (with the v dependence of 



■*/ 2 ) and suppressed for the time being), and use has 
been made of the fact that vanishes. In the propaga- 
tor Go, the initial energy Ei includes the photon energy 
w 7 (itself of order Q 2 ), i.e. E, = E 1 +E 2 +w 1 = E{ + E£, 
and the intermediate energy Ej may include, in addition 
to the kinetic energies of the intermediate nucleons, also 
the photon energy, depending on the specific time order- 
ing being considered. 

The current operators j'™' up to order n = 1, i.e. eQ, 
have been derived in Ref. In that case, the deriva- 
tion is fairly straightforward as j*- -3 -* vanishes: the low- 
est order (n = —2) contributing to j consists of the 
single-nucleon convection and spin-magnetization cur- 
rents. The situation for the charge operator is consid- 
erably more complicated, however, since n = —3 is the 
lowest order contributing to it — in momentum space, it 
is given by 



-,(- 3 ) 



(q) = e e N ,i (2ir) 3 5( Pl + q - pi) + 1 ^ 2 , (36) 



where ej\r,i = (1 + r, jZ )/2 is the proton projection opera- 
tor, q is the momentum carried by the external field, and 
the counting e Q~ 3 follows from the product of a factor 
eQ° associated with the 7AA vertex, and a factor Q~ 3 
due to the momentum-conserving 5-function implicit in 
a disconnected term of this type, see panel (a) in Fig. [5] 
below. Therefore, the operators and p^\ obtained 
from Eqs. ([34]) and ([35]) . depend on the off-the-energy- 
shell extensions adopted for and v^- 3 \ In particular, 
it appears that all of these extensions lead to a op- 
erator which i) is free of divergencies, as required by the 
absence of countcrterms at this order (eQ), and ii) sat- 
isfies /0^(q = 0) = 0. This last condition follows from 
charge conservation, 



P(q = 0) = J dxp(x) 



= P ( ~ 3) (q = o) 



(37) 



implying p("^~ 2 )(q = 0) = 0. In Sec. [V] and App. E 
we show explicitly that the off-the-energy-shell prescrip- 
tions adopted for v^ 2 \v) and v^ 3 \u) corresponding to 
v = 0, 1 do ensure that pW(j/;q) obeys requirements i) 
and ii). Indeed, we also show that the unitary equiv- 
alence extends to the OPE p(°^(^;q) — a fact already 
known [13, E3| — and TPE p^(u;q) charge operators. 



IV. CHARGE OPERATORS UP TO N3LO 

The LO contribution to the two-nucleon charge oper- 
ator in panel (a) of Fig. [51 resulting from the first term 
of the 7 A interaction Hamiltonian in Eq. ([2]), has al- 
ready been given in Eq. ([36]) . There are no NLO (e Q~ 2 ) 
contributions, whereas at N2LO there is i) a relativistic 
correction of order (Q/itln) 2 to the LO charge operator, 
panel (b) , given by 



,(-!) - 



(2/ijv,i - ejv,i) 



'A' 



x (g 2 +2iq-(Ti xKi) + 1^2 



(38) 
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FIG. 2: Diagrams illustrating one- and two-body charge oper- 
ators entering at LO (e Q~ 3 ), panel (a), N2LO (e Q _1 ), panels 
(b), (c) and (d), and N3LO (eQ°), panels (e). There are no 
NLO contributions. Nucleons, pions, and photons are denoted 
by solid, dashed, and wavy lines, respectively. The square in 
panel (b) represents the (Q/mjv) 2 , or (v/c) 2 , relativistic cor- 
rection to the LO one-body charge operator, whereas the solid 
circle in panel (e) is associated with a jttN charge coupling 
of order eQ (see text). Only one among the possible time 
orderings is shown in panels (c), (d), and (e). 



ii) a pion-in- flight term, panel (c), which, however, turns 
out to vanish when the contributions of the six time- 
ordered diagrams, evaluated in the static limit, are 
summed up, iii) a one-pion-exchange (OPE) contribu- 
tion, panel (d), which vanishes due to an exact cancel- 
lation between static irreducible and recoil corrected re- 
ducible amplitudes @. In Eq. (|38[) and what follows, q 
denotes the momentum carried by the external field, and 
ki and K; are defined as 



k, = ( P ; + Pi )/2 



(39) 



where p.j and are the initial and final momenta of 
nucleon i. Hereafter, momentum conserving <5-functions 
(q = k,-) in p(~ 3 ' and p'" 1 ' (and ki + k2 = q in the 
following expressions of two-body charge operators) will 
be dropped for brevity. We note that the power count- 
ing is different for the electromagnetic current operator, 
for which the LO term is of order eQ~ 2 (in the two- 
nucleon system), i.e. it is suppressed by an extra power 
of Q relative to p' -3 -*, and where there are NLO (eQ" 1 ) 
corrections involving seagull and in-flight contributions 
associated with OPE, which have no counterpart in the 
present case. 

The N3LO contribution illustrated in panel (e) of Fig. [2] 
is associated with the jirN coupling of order e Q origi- 
nating from the first term in Eq. ((4]) ; it gives rise to the 
vertex 



2mAfF 7r 



(40) 



for absorption (or emission) of a pion of momentum k, 
energy uif~, and isospin component a, where (2wfc) -1 / 2 is 
the normalization factor entering the normal modes ex- 
pansion of the pion field. The two-body charge operator 
follows easily by evaluating (in the static limit) the con- 
tributions of the two time-ordered diagrams: 

(o) e 9a i , n ' q ^ • k 2 

Pc ' = n ~£ ( X l ' T 2 + T 2z ) 2 + 1^2. 



2m N F* 

(41) 

In the present xEFT context, pi " 1 was derived first by 
Phillips in 2003 0- However, it is worthwhile to point 



out that the presence of an operator of the form given 
in Eq. ([41]) has been known for some time — see the 1989 
review paper by Riska [l6[ and references therein. It was 
obtained by considering the low-energy limit of the rel- 
ativistic Born diagrams associated with the virtual-pion 
photoproduction amplitude. Subsequently, calculations 
based on realistic wave functions for the A = 2-4 nu- 
clei showed that this operator plays an important role 
in yielding predictions for the A structure function and 
tensor polarization of the deuteron [l7j , and charge form 
factors of the trinucleons and a-particle (l8j , that are in 
excellent agreement with the experimental data at low 
and moderate values of the momentum transfer (q < 1 
GeV/c). These calculations also showed that the con- 
tributions due to pe*^ are typically an order of magni- 
tude larger than those generated by the Darwin-Foldy 
and spin-orbit relativistic corrections — i.e., the operator 
p(~^ above — or by vector-meson exchanges. 

There are also N3LO contributions originating from 
non-static contributions in diagrams of type (c) and (d), 
resulting from expanding the energy denominators in- 
volving pions as in Eq. ([8]). We obtain: 



m N 

x (ki-Ki-ka-Ka 



cti ■ ki er 2 • k 2 



(42) 



The contributions from diagrams of type (d) depend 
on the off-the-energy-shell prescription adopted for 
v^\v) [ijj [H El- For v = and v = 1, we find by 
direct evaluation of the relevant diagrams: 



p2V = o)=- 



g\ ai • k 2 er 2 • k 2 



4mjv Fl uj* 2 
(Ti ■ T 2 + T2,z) q ■ k 2 + 2 i (n X T 2 ) 



xk 2 -(K!+K 2 ^ 



x k 2 • K 2 + 1 ^ 2 , 



9a 



1 ^ 2, (43) 
f i • k 2 <r 2 • k 2 



OJ, 



(44) 



and it is easily seen that they are related to each other 
by the unitary transformation U^°\v), that is 



pJV) = #V = °) + 



(0), 



,(-3) 



Pd (" = 0) + * e e NA Ut°\v; p' - q/2, p) 



[/(°'^;p',P + q/2)e W4 



1^2, 



(45) 



and p and p' are the initial and final relative momenta. 
We observe that p^ + p^iv) = in the limit q = 0, as 
required by charge conservation. We also point out that 
the isovector term proportional to (ti x t 2 ) 2 in Eq. (f4"3"l) 
vanishes in the Breit frame, where pi + p 2 = —q/2 and 

p'i + p' 2 = q/ 2 - 
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V. CHARGE OPERATOR AT N4LO 

First, we note that there are non-static N4LO contri- 
butions from diagrams of type (c)-(e) in Fig. [2] — though, 
those relative to panel (e) cancel out, when the two time 
orderings are taken into account. There is also a non- 
static contribution of order (Q/mjv) 4 to the LO p^ 3 ^ 
charge operator, which we ignore in the present section. 
Here we only deal with static N4LO corrections from one- 
loop diagrams of the type represented in Fig. [3J since 
those induced by the second term in the interaction H^^n 
(proportional to the time derivative of the pion field) van- 
ish. Thus up to N4LO included, there are no unknown 
low-energy constants entering the electromagnetic charge 
operator. 

The pion-in-flight contributions illustrated in panels 
(a) and (b) of Fig. [3J involve irreducible diagrams only, 
and they are obtained by direct evaluation, in the static 
limit, of the corresponding amplitudes. We find that the 
"football" contribution — panel (a) — vanishes, while the 
"triangle" pion-in-flight operator — panel (b) — reads 



25i 
Ft 



T2,z 



qi • q 2 



(46) 



where the and w, = (q 2 + m 2 ) 1 / 2 denote the momenta 
(with the flow as indicated in the figure) and energies of 
the exchanged pions, and the integration is on any one 
of the q,'s, the remaining q/s with j ^ i being fixed 
by momentum-conserving ^-functions — as noted in the 
previous section, an overall (27r) 3 5(ki + k2 — q) has been 
dropped. 

Diagrams illustrated in panels (c)-(j) of Fig. [3J have 
both reducible and irreducible pieces. As discussed in 
Scc. lIIH the evaluation of the amplitudes is carried out re- 
taining recoil corrections to the reducible diagrams (up to 
N4LO accuracy, in this particular instance), along with 
the static (N4LO) irreducible contributions. We find that 
recoil corrected reducible contributions partially cancel 
static irreducible terms at the same order. This is dis- 
cussed in considerable detail in App. [B] Consequently, 



(/) 



2 

G?) 



Of 



th) 



(i) 



FIG. 3: Diagrams illustrating one-loop charge operators en- 
tering at N 4 LO (e Q), notation is as in Fig. [2] Only one among 
the possible time orderings is shown for each contribution. 

at N4LO the charge operator associated with diagrams 



of type (c) shown in Fig. [3J reads 

2 5 i r i 



F 4 



2(ti,,+tv) 



J l ^2 ^3 

x (Q2 • qi q2 • q3 - <Ti • q 2 X q x <r 2 • q 3 X q 2 ) 
- (ti x t 2 ) z (qi • q 2 er 2 • q 3 x q 2 

+ q 2 • q 3 <ti • q 2 x q^ , (47) 

while that arising from contributions of type (d) diagrams 
vanishes, since the integrand is an odd function of the 
loop momentum qi. For type (e) diagrams, we find: 



M 

' F 2 



(n, 



T2, 



1 



C s qi • q 2 



+ C T (tTi • qi cr 2 • q 2 + o"! • q 2 cr 2 • q! 



qi ■ q 2 (T\ ■ <t 2 ] 



(48) 



The operator vanishes due to an exact cancellation 
between the static irreducible and recoil-corrected re- 
ducible amplitudes associated with the diagrams illus- 
trated in panel (f). For type (g)-(j) diagrams, we find: 



(i) 2g A 
Pi. - ~ e -f^t T 2.. 



(£\u = o) = 



Ft 
Ft 



qi • q 2 



i 



2. 



(49) 



T2,z (qi • q 2 ) 



J \ ^2 

n, x cr 2 • q 2 x qi <T! • q 2 x qi 

,2 , ,2 

Ft 

x o"! • q 2 x qi 



At \ f W l ~ W 2 

| (Ti X 72 z / qi • q 2 



1 



2 , 



(50) 



ti\v = l) = p[l\v = 0) + 



eM^r {ti x n) 



2 Ft 

[<Ti ■ q 2 x q : <r 2 ■ q 2 x qx 



(qi ■ q 2 ) 2 ] 



LJ7 LO. 



qi ■ q 2 



1 ^2 



x (<ri + <x 2 ) • q 2 x q! 

p[ L) (v = 0) = pW(v = l)=0, 



+ 1 



M 

3 F% 
V 2 i 



t i,z (3 Cs — Gtctx- c 2 ) 



1 



(51) 
(52) 

(53) 



and a fairly detailed overview of their derivation is in 
App.E 

A few comments are now in order. Firstly, the loop 
integrals entering the expressions above are ultra-violet 
divergent. However, the total charge operator at N4LO 
is finite, since the divergencies associated with contribu- 
tions (b) and (g), (c) and (h), and (e) and (j) cancel out. 
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This is in line with the fact that there are no countcrt- 
erms at this order. In particular, we observe that the 
constraint imposed by charge conservation is satisfied, 



(!) i (!) n J (!) 

since p h + Pg =0 and pi ' 



/5- = in the limit q = 



(or ki = — IC2), while the contribution associated with 
diagram (c) in Fig. [3] can be written as (since qi = — q3 
at q = 0) 



pW(q = 0) =e 



f 



2r 2 , z (q 2 • qi) S 



+ 2 ri >2 <7i • q 2 x q x cr 2 • q 2 x q x 
- (ti x t 2 ) 2 . q 2 ■ qi <7i • q 2 X qi 
x (2 7 r) 3 (5(q 1 -q 2 -ki) + l ^2 . (54) 

ft is then seen that, in this limit, the expression above 
is opposite in sign to that of diagram (h) in Eq. (|50| for 
v = 0. For v = 1, in Eq. (|5T|) the extra terms propor- 
tional to (ti x r 2 ) z vanish by themselves at q = 0. For 
completeness, we list the configuration-space representa- 
tion of these operators in App. O 

Secondly, the charge operators p£' for v = 0, 1 arc 
related to each other by the unitary transformation U, 
i.e. a relation similar to Eq. (|45[) holds with U^ \v) being 
replaced by U^'(u), defined in Eq. (|28|) . This is easily 
verified by expressing U^\v) as 



if/W^p'-p) 



4Fi 



(3/2-n -r 2 ) 



W 2 . + uj'i 
4 4 



[(p'-p) 2 -s 2 ] 2 -4[ < T 1 .(p'-p)xs] 



1*2 ■ (P'-P) X 8] 1 - 2i [(P' - P) 2 - ^ 



X (<Ti 



where 



0-2) ■ (p' - p) x s 



w± = \/(p' - p ± s) 2 + 4 m 2 



(55) 



(56) 



The commutator [p'~ 3 - > , t/^^z/)] is seen to be identical 
to the (ti x r 2 ) 2 term on the right-hand-side of Eq. (|5T|) . 
when the pion momenta qi i2 are expressed as qi j2 = 
q/2 + p — p' ± s and s is the loop momentum. 

Thirdly, we compared the operators given above with 
those derived b y R olling et al. @ in TOPT with the 
Okubo method [19|, in order to decouple, in the Hilbcrt 
space of pions and nucleons, the states consisting of nu- 
cleons only from those including, in addition, pions. We 
find that the expressions for operators (a), (b), (c), (g), 
and (h, v = 0) are identical to those reported in Ref. 
the terms involving contact interactions in panels (d), (e), 
(i, u), and (j) were not considered by the authors of that 

paper. We should note that in p^\v — 0) the additional 
isovector piece, i.e. the term multiplied by {t\ x t 2 ) z in 



Eq. (|50|) . is missing in Ref. However, evaluation of 
the loop integral shows that it vanishes. Indeed, consider 

q 1 .(k 2 -q 1 )<x 1 -(k 2 x_q^ = j\ x2x{l _ 2x) 



qi (q 2 +m 2 r )[(k 2 -qi) +TO 2 ]- 
qi • k 2 o-i ■ (k 2 x qi) 

1 [ql + m 2 + kl x (1 - x)] 



(57) 



after making use of Feynman's paramctrization, and 
shifting the integration variables as qi — xk 2 — > qi. 
Thus, the type (h) charge operator derived in Ref. Q 
corresponds to the v = off-the-energy-shell extension. 
On the other hand, the framework used by these authors 
leads to vanishing non-static corrections to the OPE po- 
tential [13] (see also the discussion by Phillips Q in 
connection to this issue), which would imply the choice 
v = 1/2 in Eq. This suggests that pion retardation 

effects may not have been treated consistently in Ref. [H . 
We conclude by observing that for clarity's sake we have 
kept the (vanishing) isovector terms in p^\v), Eqs. (|5"0|) 
and (l5il. 



VI. CONCLUSIONS 

We have presented a fairly systematic derivation of the 
two-nucleon charge operators up to one loop (or N4LO) 
in xEFT, based on TOPT with a careful treatment of 
the non-iterative contributions extracted from reducible 
diagrams. The specific form of the N3LO and N4LO 
charge operators depends on the off-the-energy-shell pre- 
scriptions adopted for the non-static pieces in the OPE 
and TPE potentials. This ambiguity is of no import, 
however, since these OPE and TPE (non-static) poten- 
tials and accompanying charge operators are related to 
each other by a unitary transformation. Thus, provided 
a consistent set is adopted, predictions for physical ob- 
scrvablcs, such as the few-nuclcon charge form factors, 
will remain unaffected by the non-uniqueness associated 
with off-the-energy-shell effects. 

However, it is important to stress that in the present 
work we have only examined those off-the-energy-shell ef- 
fects relating to pion retardation [l(| [Uj], which arise, in 
TOPT amplitudes, from energy denominators containing 
pion (in addition to nuclcon kinetic) energies. There are, 
of course, additional non-static corrections originating 
from the non-relativistic reduction of interaction vertices 
(generated by fully relativistic Lagrangians) . Corrections 
of this type in the OPE sector for both potentials and 
charge operators have been studied in Refs. [ill EH, Hi- 
It would be interesting to extend those considerations to 
the TPE sector, and also explore the constraints, in the 
present %EFT setting, that relativistic covariance and 
power counting impose on these non-static terms of the 
potentials and electromagnetic charge and current oper- 
ators. As a matter of fact, a study along these lines, but 
dealing only with the two-nucleon potential, is that of 
Ref. ml 
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Finally, we note that the charge operators up to N4LO 
included contain no unknown low-energy constants. The 
N4LO corrections are purely isovector, and will not con- 
tribute to isoscalar observables, such as the A structure 
function and tensor polarization of the deuteron or charge 
form factor of 4 He. They will produce, presumably tiny, 
contributions to the isovector combination of the trinu- 
cleon radii and charge form factors. A quantitative anal- 
ysis of all these effects is in progress. 
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where the fields U and u M , and the covariant derivatives 
D^N and D^U, are given by 



U+U* = 



F~ = 



= T • d u TT 



2e 



(t x tt) 2 



dp+ie e N A^ + — r ■ (it x <9 M 7r) 



— A^ [tt x ( T x tv)] z + . 

— — T • O u TT -TV ■ O u 7T 



N 



2ie 
2e 



Ah{t x tt). 



— (r x tt) 2 



(A5) 
,(A6) 

(A7) 

(A8) 
(A9) 



Appendix A: Chiral Lagrangians 

In the heavy-baryon formalism [l2|, [l3j , the chiral La- 
grangians describing the interactions among nucleons, pi- 
ons, and photons are written as 



gA, F v (Ftt = 186 MeV), and e are, respectively, 
the nucleon axial coupling constant, pion decay am- 
plitude, and proton electric charge, d&, dy, d 2 $, d 2 i, 
and d 22 are (unknown) low-energy constants (LEC's), 



denotes the anticommutator, and 



and 



5 M = (i/2) 75 <7 A " / v v are the nucleon's four- velocity and 
spin operator, which in its rest frame reduce to = 
(1,0) and = (0,er/2). We have also defined 



4it = NHtv^ + g A S^)N , 
1 



c 



(3) 
■N 



2 niN 



-Art 



igA S^v v [Dft , u v ] + 



i ed 20 S tl v u F llv [t z , u p ] v p 



e y2 di + d§ r ; 

2 Mjv — c-n 



8m 2 N 



(Al) 



(A2) 



+ i e d 21 S^F^ [t z , u v \ + e d 22 [D v , F~] 



Am 2 N 



ie^°P v a S F fi(7 v°D v 



h.c. 



Fl 



N 



(A3) 



45 = -^tr[D»UD^ + ml(U + tf)] , (A4) 



e N = (1+t z )/2 , kn = (ks+Kvt z )/2 , n N = e N + n N , 

(A10) 

where k$ and Ky a re the isoscalar and isovector combina- 
tions of the anomalous magnetic moments of the proton 
and neutron, which are related to the LEC's cq and C7 
used in Ref. [l3| as k$ = c§ + 2c7 and Ky = cq. Note 
that only electromagnetic interaction terms have been in- 

( 3) 

eluded in C^^, and that the terms proportional to the 
LEC's dg and di represent corrections arising from the 
nucleon substructure (i.e., an electromagnetic form fac- 
tor). They are not relevant to our discussion here, and 
will be ignored hereafter, together with the effects due to 
the pion cloud of the nucleon. In the expressions above, 
the electromagnetic vector and tensor fields are denoted 
by A^ and P" , the isospin doublet of (nonrelativistic) 
nucleon fields by N, and the isospin triplet of pion fields 
by tt. In terms of these, the Lagrangians are expressed 
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c 



(1) 

rN 



id i 



^-T a (T ■ Vn a - -^T ■ (tt x d°Tv) 



ee N A° + —2 A [tt x (r x tt)] 2 



F 



the reducible time ordcrings, after the energy-dependent 
terms representing iterations in the Lippmann-Schwinger 
(LS) equation have been properly identified and removed, 
i.e. cancelled by terms on the right- hand-side of Eq. f)35|) . 
N (All) ^ ne kftter will be referred to below as "recoil-corrected" 
reducible contributions. As mentioned in Sec. IIIH in or- 



A 2 'j - 

7rN — 



1 



2mjv 



-TVt 



e 3A / . N 

— — [T ■ 7T + TT Z ) a 



N , 



£ (3) - 7V t 

■kN — JV 



e (2 d 2 o + 2 d 2 i - ^22 



(A12) 

<7 • (VA°) (T X a°7r), 



e(2[i N - e N ) 



[(V 2 A°)+<t-(VA°) x ^ 



<r • ^ x (VA )] 



A 



/•(2) 



1 



(9 U 7T • <9 M 7T 



2 



e ^ (tt x 9 M 7r) 



(A13) 



(A14) 



where the term proportional to eg^jF^ in the second 
line of Eq. (|A12|) is obtained Q i) by expanding the an- 
ticommutator in Eq. (|A2p as 



■ 9 A 

2 m^v 



S^if [Dp , u„], N = i 



9 A 



m N F v 



N 



. (A15) 



ii) by removing the derivative d v acting on the pion field 
via partial integration, and iii) by using the (lowest order) 
equation of motion for the nucleon field, i.e. 



' d v N 



-ee N v u A v N + 



(A16) 



to re-express the terms, which result from ii) and involve 
iv v d v N and its adjoint. In Eqs. (|Allj) - (|A14|) . we have 
retained only linear terms in the vector potential, and 
only contributions relevant for the derivation of the two- 
body charge operator up to order e Q. Application of 
the standard rules of canonical quantization leads to the 
interaction Hamiltonians listed in Sec. [Ill 
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FIG. 4: Diagrams illustrating the T (3) NN amplitudes. The 
kinetic energies of intermediate nucleons are as given. Only 
one among the possible time orderings is shown. Notation is 
as in Fig. [2] 

der to isolate these recoil-corrected pieces from those em- 
bedded into the iterated solution of the LS equation, it is 
necessary to identify the formal expressions of the N3LO 
contributions to the NN potential. The latter are ob- 
tained by retaining terms beyond the leading one in the 
expansion of the energy denominators including pion en- 
ergies, which enter both the reducible and irreducible am- 
plitudes. Note that the N3LO contributions from higher 
order chiral Lagrangians are of no interest here. Next 
section is devoted to the derivation of the N3LO NN po- 
tential, while the last two sections deal with the deriva- 
tion of the N3LO one-pion-exchange (OPE) and N4LO 
two-pion-exchange (TPE) charge operators. 



1. NN potential at N3LO 



Appendix B: Derivation of the N4LO Charge 
Operator 

In this appendix, we derive the static N4LO correc- 
tions at one loop to the electromagnetic charge opera- 
tor which follow from Eq. ([331) . The derivation of the 
operators associated with the irreducible contributions 
illustrated by panels (a) and (b) in Fig. [3] is straightfor- 
ward. However, the analysis of the reducible diagrams of 
type (c)-(j) in the same figure is more delicate, since the 
corresponding amplitudes contain (static) contributions 
originating from two different sources: one arising from 
the sub-class of irreducible time orderings for each of the 
diagrams (c)-(j), and one consisting of the "left-over" in 



One-loop contributions to the JVjV potential consid- 
ered in this appendix are shown in Fig. [4] We discuss in 
depth the results obtained for the box diagrams shown 
in panel (a) of this figure. The remaining corrections can 
be easily derived following the steps outlined here, and 
for them we only provide a listing of their expressions. 

The classes of diagrams contributing to the N3LO box 
amplitude are illustrated in Fig. The type (a) re- 
ducible, and type (c) and (d) irreducible diagrams, are 
evaluated by keeping next-to- leading (or Q°) terms in 
the expansions of the energy denominators which include 
pion energies — see Eq. ([8]). The amplitude correspond- 
ing to the reducible type (b) diagrams is obtained by 
retaining terms of order Q in these energy denominators, 
namely one order higher than for diagrams of type (a), 
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(c), and (d). 

We write the N3L0 amplitude associated with the box 
diagrams as the sum of reducible and irreducible contri- 
butions, i.e. 



r (3) ( 



,(3) 



(Bl) 



The reducible amplitude consists of LS terms plus a term 
contributing to the definition of the NN potential at 
N3LO. The latter is affected by the choice of the off- 
the-energy-shell prescription adopted for the N2L0 OPE 

(2) 

potential v\ (y). As an example, we discuss the results 
obtained with v = and v = 1. In particular, for v = 0, 
we find that the box reducible amplitude is given by 



1 



w 2 Ei-Ex- E 2 
V 3 V 4 



ViV 2 



{E l -E l ) I + {E 2 ~E 2 y 



(£1 - E[f + (£ 2 - Eft 



1 



V 2 V 3 



Ei-Ex- E 2 wi 

Ei — E\ — E 2 



Vi v 2 v 3 v 4 



,2 , ,2 



(B2) 



where E{ = Ei + E 2 = E[ + E' 2 is the initial energy of the 
system, Ej and (wj) denote, respectively, the energies 
of the intermediate nucleons and momenta (energies) of 
the exchanged pions as indicated in panels (a) and (b) 
of Fig. [51 and an integral over an unconstrained pion 
momentum is understood. In the equation above, and 
through the remainder of this appendix, we denote with 
Vi the vertices entering the diagrams. These vertices are 
implied by the interaction Hamiltonians listed in Sec. [TXJ 
For example, the V\ vertex shown in panel (a) of Fig. [5] 
is associated with the H n N Hamiltonian, and reads 



Vi 



9a <*\ ■ q 2 

-l— ^=Tl 6 



(B3) 



where b specifies the isospin component of the pion. 

The last term in Eq. (|B2|) is the N3LO recoil-corrected 
reducible contribution to the NN potential mentioned 
earlier, corresponding to the prescription v = 0. After 
resolving the spin-isospin structure implied by the ver- 
tices, one can easily recognize that the first two terms 
in ; rod represent iterations of the LS equation with 

the static, and N2LO, vi 2) \v = 0), OPE potentials 
defined in Eq. (fT7|) and (fT6|) . respectively, namely 



T. 



(3) 



27T. rod 



[y = 0) = 4°) Go („ = 0) + [y = 0) G v$ 



(0) 



VI V 2 V 3 V A 



Ei — Ei — E 2 



(B4) 



where, for brevity, the dependence upon nucleon energies 
and pion momenta is not explicitly indicated. It can 
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FIG. 5: Diagrams illustrating the recoil-corrected reducible, 
panels (a) and (b), and irreducible, panels (c) and (d), ampli- 
tudes contributing to the AW potential at N3LO. Pion lines 
with crossed (full) circle indicate that only the next-to-leading 
Q (next-to-next-to- leading Q 1 ) term in the expansion of en- 
ergy denominators, Eq. @, are retained in the corresponding 
amplitudes. See text for explanation. In panels (b), (c), and 
(d), the crossed or full circles can be either on pion one (as 
shown in the figure) or on pion two. Only one among the 
possible time orderings is shown. Notation is as in Fig. 



be inferred from Fig. [5J If the prescription 



1 is 



(2) 

considered for v\ [y] , the box reducible amplitude at 
N3LO reads instead 



T, 



(3) 



Ay = 1) = t'[ 0) Go v^{y = 1) + v%\v = 1) G 



(0) 



+ Vi v 2 v 3 v 4 

Ei -Ei- 



(B5) 



(2} 

with t4 ' (y = 1) as given in Eq. (jl8p . provided the rele- 
vant nucleon energies are considered. 

To complete the evaluation of the box amplitude, we 
need the expression of the irreducible contribution, which 
is given by 



T, 



(3) 



27r, irr 



Vi v 2 v 3 v 4 



Ei — Ei — En 



2 , ,2 



wr uj 



1 ^2 



-2 Vi Vi V 3 V 2 



Ei — Ei + E 2 — E 2 



L0fuJ 2 



E{ - Ei 



E 2 — E 2 



LOl Ll>2 



(B6) 



where, referring to Fig. [SI the first term is from the (irre- 
ducible) direct diagrams of the type shown in panel (c), 
while the last two are generated by the diagrams of the 
type shown in panel (d). The sum of the reducible and 
irreducible pieces can then be written as 

T$ {v) = 4°) Go {v) + {v) Go + «£? M , 

(B7) 

where v^J (u) for v = 0,1 are as given in Eqs. (fT^)) 
and ([20]). 

An analysis similar to that outlined above leads to the 
following expressions for the "triangle" TPE, panel (b) of 
Fig- HI an d OPE contact, panels (c) and (d), amplitudes 
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1 



(6) 



#1 



(c) 



(d) 



(e) 



(/) 



FIG. 6: Diagrams illustrating recoil-corrected OPE amplitudes contributing to the charge operator at N3LO. Only one among 
the possible time orderings is shown. Notation is as in Figs. [5] and [S] 



and corresponding potentials: 



T (3) _„(3) 



9 2 a f 1 / 
■ti -T2 / -3— f(qi ■ q2 



F 4 



! o-i • q 2 x qx)^ + E{ — 2 Ei) + 1 ^ 2, (B8) 



= ^Go4 2) M +4 2) MGo 

+ "CT,cM > 



(0) 



T {3) _ (3) _ 3^ f <ri ■ qi ^ <7i ■ qi 

i CT,d-«CT,d- ^ J -4 

x(E 1 -E 1 +E' 1 -E' 1 ) + 1^2 , 



where 



(0) r< , n 
V CT = Og + Ut (Tl ■ tT 2 , 



(B9) 



(BIO) 



(Bll) 



is the contact potential at LO, while the N3LO potential 
arising from the diagrams of panel (c) with the choices 



(2) 

v = 0, 1 for the OPE potential v), (v), is given by 



,(0) 



9a I CT i ' w ct • qi 

X (E 2 -E 2 +E' 1 -E' 1 ) + 1^2 , (B12) 
/■ (o) (0) 

4 3 > > = 1) = „g> > = 0) + / ^^(^ - id - £ 2 ) 
2^(15, -2i* -2$), (B13) 

and the nuclcon energies and pion momenta are defined 
in Fig. |U 



convenience, is obtained as 



(0) = _ 4 Vl V2 VjTT 

U>1 UJ 2 



{E 1 ^E[-E 2 +Ei) 



(B14) 



The vertex Vy W is proportional to that associated with 
the interaction Hamiltonian H J7r in Sec. [HI 



1 e e abz 



(B15) 



where the pion energies are as indicated in panel (a) of 
Fig. El It is convenient to factor out the energy numerator 
(oji — 0J2), and define V 1T: as 

V-r* = -ie e abz . (B16) 

V4 uji 0J2 

Next, we consider the amplitude associated with the di- 
agrams shown in panels (c)-(f) of Fig. [6) 

flW = [vM{v) Go P ™ Go v^M 

+P d %) . (B17) 

where p(~ 3 ' is the LO charge operator given in Eq. (|36jl . 
while p^ ^^) is the N3LO OPE contribution defined in 
Eqs. (|i§ |) -([3S p for v = 0, 1. The latter is written as 

p d °V) = $k") + p^M + 1 - 2 > ( Blg ) 

where p^jL comes from the diagrams shown in panels (c) 



and (d) of Fig. [6l while p Nl is associated with those of 
panels (e) and (f). For v = 0, we find 

Pi> = 0) = (E{ -E 1+ E 2 - Ei) , (B19) 



V Nl V x V 2 



(E[ -E 1 +E 2 - Ei) ,(B20) 



2. OPE Charge Operators at N3LO 

Before turning our attention to the N4LO corrections, 
we outline the derivation of the N3LO OPE charge op- 
erators whose expressions have been given in Eqs. (|42| 
flU]) . As discussed in Sec. IIV1 these operators vanish in 
the static limit (that is, at N2LO), while at N3LO they 
are given by amplitudes associated with the diagrams of 
the type illustrated in Fig. [5J In particular, the charge 
operator in Eq. (|4"2"j) . which we denote for later 



while for c = lwc obtain 



P<> = 1) 



2 Yl^lN { e 2 _ E i) , (B21) 



, ,3 



= = -2^^(E 2 -Ei). (B22) 

The energies E\ and E l are as indicated in panels (c) 
and (e) of Fig. O respectively, and V 1 n — eejv,i is the 
vertex implied by the interaction Hamiltonian H 7 n at 
LO. 
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FIG. 7: Diagrams illustrating the static, panels (d), (e), (i), 
and (j), and recoil-corrected, remaining panels, diagrams as- 
sociated with the N4L0 contribution shown in panel (c) of 
Fig- El Only one among the possible time orderings is shown. 
Notation is as in Figs. [2] and [5] 



3. N4LO Charge Operators 

We can now proceed to sketch the derivation of the 
charge operators illustrated in panels (c)-(j) of Fig. [3] 
The one-loop corrections of panels (c)-(e) involve a jtttt 
electromagnetic vertex, while a jNN interaction enters 
those of panels (f)-(j). We give details on the derivation 

of the operators pc and as representatives of these 
classes of diagrams. 

First, we consider the pion-in-flight term. The dia- 
grams contributing at N4LO are shown in Fig.0 As dia- 
grammatically shown in the figure, the irreducible contri- 
butions, panels (d), (e), (i), and (j), are evaluated in the 
static limit, while next-to- leading order terms in the ex- 
pansion of energy denominators are retained in the eval- 
uation of the reducible contributions. We write the total 
amplitude as a sum of these, that is 



7c 



T (i) , T (i) 

7c, red ' 7c, irr 



(B23) 



The reducible and irreducible contributions are given by 



2! 



(i) 



7c, red 



(i) 



2 (ui - uj 3 ) 
u>x oj 2 uj 3 {loi + uj 3 ) 
2 (wi - uj 3 ) 



Wl UJ 2 UJ3 (^1 + LJ3) 

[Vi V 4 , V 3 V 2 ] V 77r , 

OJl 0J 2 W3 



[Vi V 2 , V 3 V 4 ] V 1V , (B24) 
[Vi V 2 , V 3 Vi] V 77r 

(B25) 



where [...,...] denotes a commutator and pfy is the 
N3LOOPE charge operator defined in Eq. ([Bl4)l . In 
Eq. (|B24|1 . the terms multiplied by the spin-isospin com- 
binations V\V 2 V 3 Vi Vyn and V 3 ViV\V 2 V 11X come from 
the diagrams shown in panels (a)-(c) and (f)-(h), respec- 
tively, whereas in Eq. (|B25[) the first term results from 
the evaluation of the diagrams shown in panels (d) and 
(i) of Fig. [7] and the second one is from those illustrated 



in panels (e) and (j). The total amplitude is then given 
by 



T« = 4 0) GoP^+p(°»Go4 0) +pW , (B26) 



where 



,(D 



ui\ uj 2 uj 3 



[Vi v 4 , v 3 V 2 ] v r 



(B27) 



which, after resolving the spin-isospin structure implied 
by the vertices, reduces to the N4LO charge operator in 
Eq. (dZJ). 

We now turn our attention to the one-loop correction 
shown in panel (h) of Fig. [3] For this contribution, we 
distinguish among three classes of diagrams depending 
on whether the photon is absorbed before pion one, class 
A, after pion one, class B, or after pion two, class C. 
These classes are represented in Fig. [8] — the vertices and 
kinetic energies of intermediate nucleons are as indicated 
in the figure. 

We start off by discussing the result obtained for the 
class A amplitude. In Fig. [5] we show the diagrams con- 
tributing at N4LO. The irreducible diagrams, panels (f) 
and (g) of this figure, are evaluated in the static limit. 
The N4LO recoil-corrected contributions associated with 
the single (double) reducible diagrams, panels (c)-(e) of 
Fig. [9] [panels (a) and (b)], are obtained by retaining Q° 
{Q 1 ) terms in the expansions of energy denominators in- 
volving pions. The N4LO amplitude is then written as 



r (i) 

7A, re> 



'7A,irr ; 



(B28) 



where 



7S rc > = 0) 



4°) Go4 2 V = 0)GoP ( - 3) 



+ v$ {v = 0) Go p<- 3 > + «(°> Go p%v = 0) 



1 



1 



_u\uj 2 bJi lo 2 (loi + UJ 2 ) 2 Ujf U 2 {ljJl +u 2 ) 

x VI V 2 V 3 V 4 V lN 

1 2 2 



Wl W2 (Wl + W2) 2 U)j W2 Wl Wj 

x Vi Vi V 3 V 2 V lN + 1 ^ 2 . 



(B29) 



The N4LO LS terms arising from the reducible diagrams 
are listed in the first two lines of the equation above, 
where ui , v^\v = 0) and v^J (v = 0) are the LO, 
N2LO and N3LO components of the NN potential given 
in Eqs. (fTT)) . (fT6| . and (p~9|) . respectively, while the p' -3 ' 
and p ®h charge operators have been defined in Eqs. ([36]) 
and (|B19p . The last two terms in Eq. (|B29|) constitute 
the N4LO recoil-corrected contribution associated with 
the reducible diagrams. In particular, the second term 
is generated by the direct reducible diagrams of panels 
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CLASS B 
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pi 
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FIG. 8: Classes of diagrams associated with the one-loop con- 
tribution illustrated in panel (h) of Fig. [3] Notation is as in 
Fig.H 



(a)-(c) and (e), while the last one is obtained from the 
contributions of type (d). 

The irreducible amplitude from diagrams in panels (f) 
and (g) of Fig. [9] reads 



T 



(i) 



7 A, irr 



l 



i 



WlW 2 (wi+W2) 2 u\ UJ 2 (0Jl + L0 2 ) 

x V 1 V 2 V 3 V 4 V yN 

2 1 



x Vi V 4 V a V 2 V lN + 1 ^ 2 , 



(B30) 



and combining Eqs. (|B29|) and (|B30j) leads to a total 
amplitude, which can be written as 

+ «[ 2) W Go ti 0) Go p<- 3 > + w&> (i/) Go 



+ 4°) Go + ^ V) + 1^2, (B31) 



where 



-g— Vi V 3 [V 2 , Vi] V lN . (B32) 

U){ UJ 2 



The derivation of the class C amplitude (see Fig. [SJ is 
analogous to that described above. We find: 

+ p^G«vV\v)G v^+p^G«v%{v) 
+ p<®( V ) Go + p W {v) + l^2, (B33) 

where pffl is the OPE charge operator defined in 
Eq. (|B20|) . and 

p£\v = 0) = 



UJ\ u> 2 



V lN V 1 V 3 [V 2 , V 4 ] . (B34) 



Class B of box diagrams at N4LO are shown in Fig. [TO] 
The reducible amplitude, associated with the diagrams in 
panels (a)-(d), is found to be 

TXc> = 0) = 4° W~ 3) Go 4 2 V = 0) 
+ ^ 2 V = 0)G p ( - 3) Gov® 
+p$ (y = 0) Go 4 0) + 4 0) G p%\ {v = 0) (B35) 

while the irreducible amplitude, corresponding to the di- 
agrams in panels (e) and (f), amounts to 

1 



T 



(i) 



7B, irr 



,2 ,.,2 



Vi V 2 V lN V 3 Va 



,.,3, .,3 



Vi Vi Vy N V 3 V 2 + 1 ^ 2 ,(B36) 



and the class B amplitude is then given by 
T$(u) = ^ Go pHO Go ««(,,) 
+ vW{v)G 0P ^G vW 
+ p%(v)G v^+v^G oP P y (u) 



(B37) 



where 



■ 2 3' ^1 V 3 [Vi , VS,] . (B38) 



Finally, the total amplitude associated with the dia- 
gram shown in panel (h) of Fig. [3] is given by the sum of 
the A, B, and C amplitudes, that is 

p ( ~ 3) G «i 0) G v^(v) + permutations 
p^G v^(u) + v^(u)G p^ 

+ ^ X) M, (B39) 
where p^ is defined as in Eq. (|B18[I . while 

p£H") = A>) + + P [ c\v) + 1^2. (B40) 

For v = 0, pg \ and p^' 1 are as given in Eqs. (|B32I) . 
(|B38I) . and (|B34|1 . respectively, and carrying out the spin- 
isospin algebra leads to Eq. ([50)) . Similarly, for v = 1 we 



find: 










1) 




-0) + 




1) 


= /£>(" 


= 0), 




1) 


= /£V 


= 0)- 



— 3 Vi v 2 v 3 v 4 v 7JV , (B4i) 

LUi UJ 2 



(B42) 



3 V lN Vi V 2 V 3 Vi , (B43) 

OJl L0 2 
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FIG. 9: Diagrams illustrating the static, panels (f) and (g), and recoil-corrected, remaining panels, class A diagrams contribut- 
ing at N4LO. In panels (b)-(e), the crossed and full circles can be either on pion one (as shown in the figure) or on pion two. 
Only one among the possible time orderings is shown. Notation is as in Figs. [2] and [5] 



from which it follows that 
Ph(v = 1) 



Ph(^ = 0) 
2 



[Vi Vi V 3 V 4 , V lN ] + 1 ^ 2 



(B44) 



and simplifying the spin-isospin structures leads to the 
operator p^\v — 1) given in Eq. ([ST]) . 

Below we list the expressions for the amplitudes as- 
sociated with the remaining N4LO one-loop corrections 
illustrated in Fig. [31 in particular, referring to panels (d), 
(e), (g), (i), and (j) of this figure (type (f) operator van- 
ishes as pointed out in Sec. [Vj we obtain: 



T. 



(i) 



(0) 
CT 



g ' 

p(~ 3 ^ G Q Uq^ G (u) + permutations 



(B45) 
(B46) 

(B47) 



CT 

+ [p^ Go v&(v) + «g> » Go p ( ~ 3) ] +1^2 



p d °\v)G v^ + v^G p d »'(v) 
}^ d G p^+p^G Q v^ d 



(°) JL .,(0) 



,( )/ 



(B48) 



(B49) 



where pe , pg , and p^' are given in Eqs. ((48]). ([39)). 

and (|53p . respectively, while p? with i/ = 0, 1 vanishes. 
In the equations above the LS terms involve the LO con- 
tact, Eq. ([BTT]) . and N2LO OPE components of the NN 
potential, Eqs. ([IB")) and (fTS)) . as well as the N3LO po- 
tential in% derived in Sec. IB II In addition, the p( -3 ' is 
defined in Eq. ([36]) . while the N3LO OPE charge opera- 
tors, p(°>, are listed in Sec. 



Appendix C: The N4LO charge operators in r-space 

The configuration-space representations of the charge 
operators at order n = — 3,...,0 are well known 1 1 81 ] - 



Here we obtain those corresponding to p( 4 > for v = 
only. Note that pi 1 ', py, and p^\v = 0) vanish. Next, 



consider 
S 4 



2 9 a 

pb (q) =e -^r r 2. 



e tk2r2 <5(k 1 +k 2 -q) 



ki,k 2 



^N<li+<l2- k i) +1^2 



-e^4 
F 4 

X 7T 

1^2 , 



r 2 , z e^ r2 [V/ OT (r)] • [V/ OT (r)] 



(CI) 



where in the second line r denotes the relative position 
r = ri — r 2 of the two nuclcons, and 



1 



1 e" 



p a + m 2 4ir r 



(C2) 



Of course, the expression above is ill-behaved in the limit 
of vanishing internucleon separations, and needs to be 
regularized. This can be accomplished by replacing 



Mr) -> Mr) 



<?a(p) 
p 2 + m 2 



(C3) 



and in applications so far j22| the cutoff function has been 
taken as C\(p) = exp(— p 4 /A 4 ). Similarly, we find: 

p^=~e 2 ti 



F 4 



2 (ti, Z + T 2 , z ) 
X [eap^CXftu <Tl tCt (7-2, A + <5/3 7 <5pi/] 
— (Tl X T 2 ) 2 eA^i/ CT2.A — e a; g 7 (5^,, CTl,a] 

^[d h pd 2 ^c^- R K(r)] [d y d u Mr)} , (C4) 
p« =-e 2 -0{r^+T 2 . z )c^8{v)I^) (C5) 

P^^U.*^"' 1 [V/.(r)] ■ [V/ W (r)] 



F 4 
1^2 , 



(C6) 



% 9 A iq-R 



F 4 



2 ti,z e Q/ 3 7 eA A ,iy o"i, Q cr 2 ,A 
+ 2 r 2:Z SpyS^ — (ti x x 2 ) z e a( g 7 <5 M1/ cti :Q 
fyfy/U*-)] [a 7 ^/ ff (r)] + 1^2, (C7) 

P? = e ^ r L* ° iq ' R + 1 - 2 > ( C8 ) 
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FIG. 10: Diagrams illustrating the static, panels (e) and (f), and recoil-corrected, remaining panels, class B diagrams con- 
tributing at N4LO. In panels (b)-(d), the crossed and full circles can be either on pion one (as shown in the figure) or on pion 
two. Only one among the possible time orderings is shown. Notation is as in Figs. [2] and [5] 



where R = (r x +r 2 )/2 denotes the two-nucleon center-of- 
mass position, the functions h„(r) and f„(r) arc defined 
as 



hJr) 



1 

8^ 



1/2 



-Lr 



dye iycir 



-1/2 



1 



1 e" 



lp {P 2 + ml) 2 8tt rriir 



where 



L= yjml + q 2 (1/4 - y 2 ) 



(C9) 



(CIO) 



the gradients (or partial derivatives) V, Vi, and V2 act 
on the variables r, ri, and r2, respectively, and 



dxe 



Cs 5 a p — CT{2&l,a 0"2„ 



cri ■ er 2 Sap) [d a f-K{x)}[dpf 7 ,{x)\ , (Cll) 



Regularized expressions are obtained via the replace- 
ments 



M')">Mr) = dye iv *' [ e ip 

7-1/2 Jp 



/2 

C a (p) 



(p 2 + L 2 ) 2 



(C12) 



ip-r CA(p) 



Mr) -> A(r) = / e l P r , 2 , (C13) 



<5(r)^. 9A (r) = / e !; P r C A (p) 



(C14) 



Lastly, we observe that i) + is proportional to 
/(q) — J(0), and this quantity remains finite for any q 
value; ii) the requirement = at q = is satisfied 
also when the cutoff A is included. 
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